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We study properties of dilute polymer solutions which are known to depend strongly on polymer 
elongation. The probability density function (PDF) of polymer end-to-end extensions R in turbulent 
flows is examined. We demonstrate that if the value of the Lyapunov exponent A is smaller than the 
inverse molecular relaxation time 1/r then the PDF has a strong peak at the equilibrium size Rq and 
a power tail at i? 3> -Ro- This confirms and extends the results of Q. There is no essential influence 
of polymers on the flow in the regime Ar<l. AtA>l/T the majority of molecules is stretched to 
the linear size Rop S> -Ro- The value of Rop can be much smaller than the maximal length of the 
molecules because of back reaction of the polymers on the flow, which suppresses velocity gradients 
thus preventing the polymers from maximal possible stretching. 

PACS numbers 83.50.Ws, 61.25.Hq, 47.27.-i, 05.40.-a 
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Dynamics of dilute polymer solutions is an important 
subject both from theoretical and practical points of 
view. Possible applications rely mainly on the fact that 
low concentrations of polymer molecules can lead to sub- 
stantial changes in hydrodynamics. The most striking 
effect related to polymers is probably the so-called drag 
reduction in turbulent flows. A consistent explanation of 
this effect is a long-standing question One believes 
that the drag reduction is related to the effective increase 
of the viscosity due to the presence of polymers |^ . Here 
we address some aspects of this phenomenon. 

An important underlying property of polymers is their 
flexibility At equilibrium, a polymer molecule 

coils up into a spongy ball with a radius Rq. For dilute 
solutions with concentrations n satisfying hRq ^ 1, the 
influence of equilibrium size molecules on hydrodynamic 
properties can be neglected. When placed in a flow, the 
molecule is deformed into an elongated structure of el- 
lipsoidal form which can be characterized by its end-to- 
end extension R. Since the number N of monomers in 
a long-chain polymer molecule is large, R can be much 
larger than Rq. It explains why minute amounts of poly- 
mers can produce an essential effect. It was shown in 
101 that in sufficiently intensive flows polymer molecules 
get strongly extended due to stretching. This is the key 
mechanism providing an essential back reaction of the 
polymer molecules on the flow. 

Here we consider turbulent dynamics of polymer solu- 
tions. We assume that R is always much smaller than 
the viscous length of the turbulent flow, r^. Therefore, 
molecules can be treated as immersed into a spatially 
smooth external velocity field Q. In this case the dy- 
namics of polymer stretching is determined only by the 
gradients of the velocity. Since the gradients in turbu- 
lent flows are correlated at the viscous length, all the 
molecules inside regions with size of the order of ry are 



subject to the same gradient, and therefore are stretched 
coherently. As long as one can neglect the hydrodynamic 
interactions between molecules, the problem is reduced 
to dynamics of a single molecule. 

We investigate the behavior of polymer molecules with 
the extensions R satisfying Rq <^ R <ti i?max, where 
Rmax is the maximal size of the polymer. Random walk 
arguments show that the entropy of such molecules is 
quadratic in R in this interval, which leads to Hooke's law 
(see e.g. Ref. [^). That is why one can expect a linear 
dynamics of the molecules. Even though hydrodynamic 
interactions of monomers make polymer's dynamics in- 
herently nonlinear, the interactions can be neglected for 
elongated molecules. This expectation is conflrmed by 
recent experiments with DNA molecules where an ex- 
ponential relaxation of a single molecule was observed. 
Numerics and theoretical arguments presented in Ref. j|] 
also show the linear character of the molecule dynamics 
for Rq R ^ Rmax- In experiments a number of the 
molecule eigen modes has been seen. We will take into 
account only the mode which has the largest relaxation 
time T, because the other modes are harder to get excited 
in turbulent flows. 

A starting point of our theory is the dynamic equation 
for the vector R which can be defined, say, via the inertia 
tensor (per mass of a monomer) RaRp of the elongated 
molecule. Then R determines the orientation and the 
largest size of the molecule. We assume the following 
dynamic equation for the vector (cf. Refs. |^,^) 



— Ka — UpW pVa , 

at T 



(1) 



where r is the relaxation time. The velocity gradient 
must be taken at the molecule position. The role of non- 
linearity in the extended equation for R (and in the sys- 
tem of equations for N coupled beads) is examined in 
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Ref. p. For our purposes this non-hnearity as well as 
the thermal noise is irrelevant (see the discussion below). 

To get rid of inessential degrees of freedom responsi- 
ble for the orientation of the molecule we write R = Rn 
passing to the absolute value R of the vector R. Then 
we obtain from Eq. (§) (cf. 0) 



dt ^ T 



dUa „ 



R = Ro exp(p) , ( = naUpViBVa 



(2) 
(3) 



We see that the evolution of p is determined by the scalar 
function C which is a functional of the velocity field. 

For turbulent flows where the velocity randomly varies 
in time one should use a statistical approach. A natural 
first step is to take the polymers being passively embed- 
ded into the fluid, disregarding their back reaction on 
the flow. We will demonstrate that there is a wide re- 
gion of applicability of this approximation. Neglecting 
the back reaction we can treat the velocity dynamics as 
independent of polymers. Then defined by Eq. (^ is 
independent of p. We will not specify the velocity statis- 
tics. Irrespective of its character one can use the large 
deviation theory (see e.g. Ref. devoted to different 
aspects of Lagrangian dynamics in turbulent flows) . The 
scheme presented below is valid for any random flow. In- 
tegrating Eq. (0) we get 



p{t) = pq + z- 



t 



z^ dt'Cit'), (4) 
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where po is the value of p at t = 0. One should keep in 
mind that the expression (^) for p is correct if one can ne- 
glect the presence of the boundaries Rq and i?max where 
Eq. (|l|) is violated. 

The integral z in Eq. (|^) possesses some universal 
properties for times much larger than the correlation time 

of the random process For turbulent flows r^; can 
be estimated as the characteristic time of the Lagrangian 
motion on the viscous scale, which coincides with the 
characteristic inverse strain on this scale. For t ^ the 
variable z can be considered as a sum of a large number 
of independent variables. Then in order to establish the 
statistics of z for fluctuations near its mean value one 
can use the central limit theorem. If we are interested in 
large deviations from the mean, a more general formula- 
tion is needed (see e.g. [^^). Namely, the PDF of z 
can be written as the homogeneous function 



: exp 
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(5) 
(6) 



"The entropy density" S' is a functional of the velocity 
statistics. It is impossible to calculate S without knowing 
the statistics explicitly. Fortunately, only general proper- 
ties of S (such as positivity and convexity) are needed for 



us. The central limit theorem is reproduced by Eq. (|^) 
if to consider a vicinity of the entropy maximum where 
S{x) w a;V(2A). The constant A defined in (|) is the 
principal Lyapunov exponent of the turbulent flow, which 
is the average logarithmic rate of growth of the distance 
between two initially close Lagrangian trajectories. 

As follows from Eq. (^, G(t,z) determines the con- 
ditional probability that p{t) has the value po + z — t/r 
provided p(0) = po- Therefore one can write the equation 



rit,p)= dpogit,p~ po + t/T)rio,po) 



(7) 



for the PDF V{t,p). In the stationary case V is t- 
independent and Eq. (|^) can be treated as a relation 
determining the PDF. Writing V as the Laplace integral 
V{p) = J da exp{—ap)V{a) we observe that the convo- 
lution in Eq. (Q) becomes a product and the equation 
can be easily resolved. The value of V{a) is non-zero if 



dx 
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(8) 



Apart from the trivial solution a = this equation de- 
fines a uniquely. 

Since i S> t^, one can use the saddle-point approxima- 
tion in calculating the integral (||). It gives the condition 

a^S'iP+l/T-\), (9) 

where /3 is the saddle-point value of the ratio (p — po)/t. 
Equating the integral in the left-hand side of Eq. (H) 
(calculated in the same approximation) to unity we get 
the equation for /? 

S (^f3~ X + ^^ - pS' (^P- X + ^^ =0. (10) 

It is important that (3 is independent of t and p. Solving 
Eq. (|l^) and substituting the result into Eq. (^) we find 
the exponent a. The trivial solution /3 = A — l/r of Eq. 
(|l0|), corresponding to a = 0, should be discarded. 

We conclude that a single component V{a) is non-zero 
and therefore V{p) oc exp(— ap). Recalculating this dis- 
tribution of p into that of R we obtain the power tail of 
the PDF of the molecule size R 



V{R) - R'^R- 



(11) 



For positive a the normalization integral J dRV{R) is 
determined by small i?, which means that the majority 
of molecules has nearly equilibrium size. On the contrary, 
the normalization integral diverges at large i? if a < 0. 
Then the majority of molecules is strongly stretched. 

Another way to obtain the result ( pT| ) is to consider 
the typical fluctuation making the largest contribution 
into the tail of the PDF. Starting from a nearly equilib- 
rium shape, that is from po ~ 1, the velocity stretches 
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the molecule up to p ^ 1. The contribution of fluctua- 
tions with stretching period t is equal to G{t,p + t/r). 
It has a sharp maximum at time determined from 
dQ{t, p + t/r)/ dt = 0. This condition gives = p/[3. 
The probability density is thus dominated by fluctua- 
tions with stretching period t^. It is proportional to 
Git*, p + t^/r) which reproduces Eq. (|l^) with a given 
by Eq. (||) . Note that the characteristic value of the ve- 
locity gradient for the relevant fluctuations is given by 
C ~ p/i* + 1/t and is of the order l/r. 

Let us establish the dependence of a on the control 
parameter, which is the strength of velocity fluctuations 
at the viscous length measured by the Lyapunov expo- 
nent A. As A tends to zero, the function S{x) contracts 
to X = and therefore a tends to infinity, which implies 
strong suppression of the tail. It is quite natural since 
in a weak flow the molecules are only weakly stretched. 
Note that even for intense flows the Lyapunov exponent A 
is suppressed in the regions where the rotation rate dom- 
inates the strain rate. As A increases, the exponent a de- 
creases and at a certain level of fluctuations approaches 
zero. If A is close to l/r then one can use the quadratic 
approximation for S which leads to the law 

/3^i-A, «=1(^-a). (12) 

We see that a changes its sign at A = l/r. Thus the 
majority of molecules becomes stretched when A > l/r. 
This can be interpreted as the criterion for the coil- 
stretch transition in turbulent flows discussed in . 

We can use Eq. (|^) only if p and po belong to the 
asymptotic region between zero and Pmax, where Eq. (|l|) 
is valid. The saddle-point approximation used above 
gives po — p — (3t. Thus the above scheme works only 
ii t < t* — pl(i (here we assume /3 > 0, i.e. a > 0). 
Then the polymer molecules spend most of the time 
fluctuating near the equilibrium shape, occasionally get- 
ting stretched by strain fluctuations which overcome the 
elastic reaction. The fluctuations leading to a given R 
have the duration k, p/ (3. Since /? tends to zero 
when A — > 1/t one should observe a critical behavior 

oc |A — 1/t|^^ in accordance with Eq. ([l^). We see 
that in the vicinity of A = l/r the time t» is much larger 
than tq which justifies our scheme. Similar considera- 
tions are valid for a < 0. 

One can generalize the scheme taking into account a 
number of molecular eigen modes. Since the critical value 
of A is determined by the inverse relaxation time, then 
in the vicinity of the critical value corresponding to the 
principal mode, the other modes are at most weakly ex- 
cited. However, they can be important at larger A. 

The rest of the paper is devoted to the discussion of the 
back influence of the polymers on the flow. A consistent 
investigation should be based on the complete system 
of equations coupling turbulence with polymers. One of 
these equations is the modified Navier-Stokes equation 



{dt -f- vW)vc, = -V„p + vW'^v,, + WfsIVo,p , (13) 

where 11^^ is the polymer contribution to the stress ten- 
sor. Equation ( |l^ ) should be supplemented with the 
equation describing dynamics of IIq,^. In the considered 
case IIq^ can be defined as a sum of stresses of polymer 
molecules in a volume divided by the mass of the fluid 
inside the volume ||l^ . We are interested in the situation 
when the molecules are strongly elongated. Then due to 
Hooke's law the stress of such molecule is proportional 
to RaRp. Next, taking the volume smaller than the vis- 
cous length we deal with coherently elongated molecules. 
Therefore the stress tensor can be written as 

IIa/3 = Ho e:>q){2p)nanp , (14) 

where is a unit vector, IIq exp(2p) is the principal 
eigenvalue of IIq,^ and the elongated molecules corre- 
spond to p > 0. Then from Eq. (^ we get the same 
Eqs. (^) for p and n, where d/dt should be understood 
as the material derivative dt + vV. Thus the velocity 
V is now coupled to p and n via Eqs. (^3|jlJ). Note 
that the constant Hq in Eq. ( p^ is proportional to the 
concentration of the polymer molecules. 

Let us consider the PDF of R not assuming that the 
flow is unperturbed by the polymers. We start from the 
case Ar < 1. One recovers Eq. ( pT| ) if the back influ- 
ence is small, i.e. 11 <C J^Vu for the relevant fluctuations 
characterized by Vw ^ l/r. Since 11 oc R^, the polymer 
contribution in the stress tensor grows with R and the 
inequality 11 <^ vSIv is violated for the molecules with 
^ 1^ ^^back- The value of i?back can be found from the 
estimate vJt ~ Ilo^back/^o- F^^' ^ ~ ^back the back 
reaction switches on and suppresses the velocity fluctu- 
ation. Hence, the probability of fluctuations producing 
R > i?back is small and hence at i? > i?back the PDF 
decays much faster than prescribed by Eq. (pi]). 

Now we study the case A > l/r. For R <C i?back the 
polymer stress tensor is small and as explained above ve- 
locity is decoupled from the elastic degrees of freedom. 
Since stretching is stronger than the elastic force, R grows 
for any typical velocity realization. On the other hand, 
at i? ^ ^back the polymer stress influences the velocity, 
suppressing it strongly for sufficiently large R. This leads 
to a decrease in R. Therefore, the majority of molecules 
has sizes near an optimal size i?op > ^back- The PDF 
is an increasing function of i? at i? < i?op and decays 
fast at i? > Rop- In this state velocity gradients can 
be estimated as l/r ||l^. This can be proven e.g. by 
averaging Eq. (|^). It means that the Lyapunov expo- 
nent of the solution is smaller than that of the solvent at 
the same energy input. The energy dissipation is related 
mainly to the polymer stress tensor and hence Rop grows 
as the input of energy increases [Q. The effective vis- 
cosity defined as the proportionality coefficient between 
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n and Vv, also grows. Note that there exists an in- 



terval where Ron Rn 



This contradicts the widely 



accepted view that at some level of turbulent fluctuations 
there is a sharp transition between the state where most 
of the molecules have i? ~ i?o and the state where all the 
molecules are stretched up to i?max- 

We conclude that for a > (i.e. A < 1/t) the end- 
to-end extensions of the majority of molecules are of the 
order of the equilibrium size Rq, and there is no essen- 
tial contribution to the stress tensor. For a < (i.e. 
A > 1/t) extensions of most of the molecules are of the 
order of i?op ^ Rq- Then the polymer stress tensor 11 
is estimated as HqRIp/Rq- Its value can be much larger 
than the viscous contribution ^A [p^ . 

The above analysis implies that i?back <C i?max since 
at R ^ Rmax one must consider non-linear corrections 
to Hooke's law and hence to Eq. (|^). The condition 
^back ^ -Rmax IS realized for sufficiently high concentra- 
tions of polymers. Then the fluid displays non-Newtonian 
properties. When most of the molecules are stretched up 
to i?max but the back reaction is not switched on, one 
has a Newtonian fluid whose properties do not differ sig- 
nificantly from the properties of the solvent. This is the 
case for very dilute solutions where 11 <^ vVv. 

Since a turbulent flow is multiscale, the real picture is 
more complicated. We have shown that the main char- 
acteristics of the flow that determines the behavior of 
polymer molecules is the Lyapunov exponent A, which 
is defined at the viscous scale. Hence, the dynamics of 
a molecule is sensitive to the fluid motion at the vis- 
cous scale whereas the velocity varies over a wide inter- 
val of scales. Therefore the Lyapunov exponent varies 
in time and space over scales from the inertial interval. 
We thus have an "intermittent picture" : In the regions 
where A < l/r one deals with a Newtonian fluid with the 
viscosity v of the solvent, whereas in the regions where 
A > 1/t the polymers are strongly stretched and the ef- 
fective viscosity can be much larger than v. 

As the Reynolds number Re increases, the relative vol- 
ume of the regions with A > l/r increases and the aver- 
aged (over space) viscosity grows. The average value of 
Rop also grows. After it has reached the value of the order 
^max, the back influence cannot grow anymore. It means 
that the effective viscosity first grows and then decreases 
back to the solvent value v. Note that the effective vis- 
cosity varies smoothly without onset. As a consequence, 
the drag reduction also varies smoothly with Re, having 
a maximum at some intermediate value. Experiments 
seem to confirm our picture (see, e.g., [0). 

To avoid misunderstanding let us stress that we con- 
sider conventional turbulent flows which have the iner- 
tial interval of scales. In the inertial interval the polymer 



back reaction is small compared to the non-linear term 
in Eq. (p^. In principle, in some region of parameters 
the back reaction can be stronger than this non-linearity 
everywhere and the properties of the fluid are drastically 
different |Q. This case requires a separate analysis. 
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